Klein tunneling in pn and npn junctions of silicene 
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We investigate Klein tunneling of pn and npn junctions made from silicene, Si analogue of 
graphene. The conductance shows the distinct gate-voltage dependences peculiar to the topological 
and non-topological phases, where the topological phase transition is caused by external electric 
field. Namely, the conductance is suppressed in the np regime when the both sides are topological, 
while in the nn regime when one side is topological and the other side is non-topological. Further- 
more, we find that the conductance is almost quantized to 0, 1 and 2. Hence we can use a silicene 
pn junction as a digital field-effect transistor. Our findings will open a new nanoelectronics based 
on silicene. 

PACS numbers: 73.23.-b, 72.80.Vp,73.40.-c 



I. INTRODUCTION 

Technology fabricating one-atom thick systems has 
been rapidly developing since the appearance of mono- 
layer honeycomb carbon (graphene) .^Recently, Si analog 
of graphene (siliceneJpMH has been synthesized. The low- 
lying excitations of the monolayer honeycomb systems 
are Dirac fermions. Due to spin-orbit interaction (SOI), 
the Dirac fermions become massive, i.e., the energy bands 
have a gap. These massive Dirac fermion systems lead 
to a quantum spin Hall (QSH) insulator, which is orig- 
inally proposed in a grapheneP^ Unfortunately, SOI of 
graphene is tiny so that the QSH effect in a graphene has 
not been experimentally observed. SOI of silicene is, in 
contrast, thousand times larger than that of grapheneP^ 
The QSH effect in silicene is experimentally accessible. 9 

Transport of Dirac fermions shows various anoma- 
lous behaviors. A massless Dirac fermion never shows 
localization^ owing to the topological termPJ] Further- 
more, Klein tunneling occurs in Dirac fermionsP^ Es- 
pecially, for the massless case, a perfec t transmission 
through any barrier potential occurs 1121141 This originates 
from the topological phase (Berry phase) 7rP^ It is dim- 
cult to observe the Klein tunneling for elementary par- 
ticles, due to a huge mass gap. In contrast, it is pos- 
sible for semiconductors in which the band gap is small 
(~ leV). Actually, the Klein tunneling has been observed 
in grapheneP^ The topologically protected transport is 
robust against perturbations. Thus, the systems sup- 
porting Dirac fermions can exhibit unique charge and 
spin transport. They have a potential to be a new elec- 
tronics and spintronics device. 

Among such Dirac fermion systems, silicene has an- 
other advantage: the band gap is controllable by ap- 
plying an external electric fielcP^ owing to the buckling 
structure.^ A bila yer gr aphene also has an electric-field- 
tunable band gappStll] Silicene, differently from a bi- 
layer graphene, shows the topological phase transition 
with tuning the band gap. If one applies an electric field 
whose energy is stronger than SOI, the topological phase 



transition from a QSH to non-topological insulators oc- 
curs. Also, silicene realizes various topological insulators 
by exchange interaction,^ photo-irradiatiorPSl and anti- 
ferromagnetic order P3 



These characteristics can be useful for device applica- 
tions. The most fundamental electronic device is a field- 
effect transistor (FET). FET made by silicene has an 
advantage that it has a large band gap due to SOI com- 
pared with graphene which is a zero gap semiconductor. 
In addition, the tunable band gap and the topological 
phase transition of silicene by external electric field may 
lead to a new feature for FETs. 



Charge transport properties of a silicene nanoribbon 
has been studied.^ Spin transport has been also studied 
in a bulk silicene junction under Zeeman fieldP^l On the 
other hand, we focus on the charge transport in the bulk 
silicene. In this paper, we analyze the transport proper- 
ties of pn and pnp junctions made of silicene. Control- 
ling the conductance by tuning the gate voltage, we find 
that i) the gate-voltage dependences of conductance in 
the topological and non-topological phases are quite dis- 
tinct, and ii) the conductance is almost quantized to 0, 1 
and 2. The former is an evidence for the existence of the 
topological phase. The latter enables us to use silicene 
as a FET with three digital values of conductance. Our 
results will open a new way to future nanoelectronics. 



This paper is composed as follows. In Sec. [TTJ we re- 
view the bulk properties of silicene. In Sec. |III| we inves- 
tigate the pn junction of silicene. First we calculate the 
one-dimensional model with and without the Rashba in- 
teraction. Next we calculate the two-dimensional model. 
The conductance of the two-dimensional model is almost 
the same but smeared compared with that of the one- 



dimensional model. In sec 
junction of silicene. Section [V 



IVJ we investigate a pnp 
is devoted to discussions. 
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II. BULK PROPERTIES 

First, we show the bulk properties of silicene. Hamil- 

tonia n of a silicene in the vicinity of the K and K' points 
readsPZESI 



H(k) = hv F (k x T x - k y T y r) z ) - XsoT z cr z i] z 
+ a\RT z (k x a y - k y <j x )i] z + £E z t z , 



(1) 



where <7j, 7$, and rji are the Pauli matrices for the spin 
(t and 1), sublattice (A and B sites) pseudospin, and 
valley (K and K' points) spaces, respectively, a denotes 
the lattice constant. Aso is the intrinsic SO coupling 
constant that triggers the topological phase transition 
from the non-topological to topological insulators (See 
below). The sublattice-dependent Rashba SO interaction 
Ar, also appears due to the buckling structure of silicene. 
And also, the mass term £E z t z shows up because of the 
buckling structure with an external electric field E z along 
the z-axis. Hereafter, we set H = 1. 

Here we review the topological phase diagram of sil- 
icene. The Dirac mass m for the K point (r\ z = 1) in the 
bulk silicene is given by 



m = -XsoT z a z + £E Z 



(2) 



The system is a topological insulator for Aso > £\E Z \, 
while it is a non-topological insulator for Aso < £\E Z \. 
Aso = £\E Z \ is the critical point, where the energy gap 
closes. The sign-change of the mass term means inver- 
sion between the conduction and valence bands. Corre- 
spondingly, direction of r and <x of the conduction and 
valence electrons changes. Namely, the sublattice and 

spin states is given by \-\ — ), | ), |++), and | — h) in 

descending energy order for the topological phase. On 
the other hand, the second and third states are inter- 
changed for the non-topological phase, i.e., |-| — ), |++), 

| ), and | — h). Here |a/3) denotes the eigenstate with 

t z = asgn(E z ) and a z — f3sgn(E z ). This energy level 
scheme is shown in Fig. [I] As one switches off E z , the 

two states with |H — ) and | ) (with |++) and | — h)) 

are degenerated (ei = e 2 in Fig- 

In addition, we show the energy dispersions for E x = 
and E z ^ in Fig. [2] The energy E±(k) in the bulk is 
obtained to be 



El(k) = uffc : 



+ ^±^/A2 + a 2A2fc 2 + IE Z \ 



(3) 



with k = (k x + ky) 1 / 2 . The corresponding eigenvector 
u±(k;E z ) is also obtained analytically (Appendix [A]). 
The energy bands are doubly degenerated for E z = 
due to the inversion [T z a z H(k)a z T z = H(—k)] and 
time-reversal [T y a x H* (k)a x T y = H(—k)\ symmetries de- 
fined within each valley. In contrast, there is no spin- 
degeneracy for E z ^ since E z breaks the time-reversal 
symmetry defined within each valley. 
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FIG. 1. Energy level scheme for the topological (TP) (£\E Z \ < 
Aso) and non-topological (NTP) | _E Z | > Aso) phases at the 
K point, ja/3) denotes the eigenstate with t z — asgn(E z ) and 
a z =/3sgn(E z ). e 1 = A SO + t\E z \, e 2 = |A SO - l\E z \\. 
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FIG. 2. Energy dispersions in the bulk silicene for E z — (a) 
and E z = 0.5Aso (b). The parameters of the system are taken 
as follows, vp/a — 1.4eV, Aso = 3.9meV, Ar = 0.7meV. 



III. SILICENE PN JUNCTION 



In this section, we investigate a charge transport in a 
silicene pn junction, which is illustrated in Fig. [3j 




charge 
current 



FIG. 3. Schematic of a pn junction. 
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A. One— dimensional case 

1. Formalism of the scattering problem 

Firstly, we consider a pn junction of one-dimensional 
silicene. The Hamiltonian is given by 

H(x) = -iv F d x r x - XsoTzVzVz + £E z 9(x)t z 

+ iaX R d x a y T z r]z + V9(x). (4) 

Hereafter, we focus only on the K point (r\ z = 1). The 
same analysis is applicable to the K' point. 

We solve the scattering problem of the pn junc- 
tion. T he calcu lation is done by employing theo ries for 
graphenePHlMl anc i the Kane-Mele model! 21122 ' j n t ^ c 
incident side (x < 0), an external electric field is not 
applied and hence the energy bands are doubly degen- 
erated. As a result, there are two incident states for a 
fixed incident energy Ep. Wave function ip±(x) of the 
scattering state with the incident energy being E-p has 
the form as 

t/ J± (-0) = u ± (fci; 0) + r±+u+(-*i; 0) + r±_ w_ (-fe i; 0), 

(5) 

^±(+Q)=t±+u+(g + ;£,) + t±_u_(g_;.E»), (6) 

where the subscript of ip±(x) denotes the spin state of 
the incident state w±(fci;0). The first term u±(ki;0) of 
Eq. ([5]) denotes the incident state. The other two terms 
correspond to the reflected states with the same (±) and 
different (spin-flip, =f) spin states. Momentum k\ of the 
incident electron is given by 



fci = sgn(£' F ) 1 



A so 



a*Xl 



(7) 



The sign of k\ is determined so that the group velocity of 
the incident state is positive. Note that the incident state 
must be a propagating mode; E F > Ag Q . Otherwise, the 
corresponding conductance is zero, by definition. On the 
other hand, momentum q± of the transmitted electron is 
obtained by solving the following equation; 



(E F - V) 



2 2 
v F q ± 



± 



cPXiqi+lEz 



x = 0. Since charge current is conserved, the following 
relation holds. 



dH 



d(-id x 



*<o d{-id x 



^±(+0), (10) 



a;>0 



where dH/[d(— id x )} is the velocity operator. The above 
relation is reduced to 



V>±(-o) = V±(+o). 



(11) 



Solving this, one obtains the reflection and transmission 
coefficients. From the reflection coefficient r a p , the trans- 
mission probability T± is given by 



7± = l-El r ±/3l 

|8=± 



(12) 



Charge conductance G in the one-dimensional system is 
defined by 



G=-(T++T_). 



(13) 



2. Charge transport asymmetry in the nn and pn regimes 

We show results on the conductance of the one- 
dimensional pn junction in Fig. [4] The horizontal axis 
is (Ep — V)/Ep, where Ep — V corresponds to the Fermi 
energy in x > measured from the charge neutrality 
point. The vertical axis is £E z /Xso- Only the region of 
£E z /Xso > is shown since the transmission probabil- 
ity is symmetric with respect to E z = (See Appendix 
[Bj. A pn junction with two doped topological insula- 
tors (TP/TP) is realized for l\E z \ < A SO - On the other 
hand, that with doped topological and non-topological 
(TP/NTP) insulators is realized for £\E Z \ > A so . Clearly 
seen from Fig. |4j there is no qualitative difference be- 
tween the conductances with [Figs. |4|a) and (b)] and 
without [Figs. [i|c) and (d)] the sublattice-dependent 
Rashba SOI A R , 

For Ar = 0, one can obtain a simple formula for the 
reflection coefficient. The reflection coefficient r a with 
a — ± being the z-component of spin of the incident 
(8) electron is given by 



Since the group velocities of the transmitted electrons 
should be positive, the following relation is satisfied for 
the propagating mode; 



sgn(g±) = sgn(£ F - V) . 



(9) 



For the evanescent mode, on the other hand, Img-t > is 
satisfied. And note that when \Ep — V\ < |Aso — £ \Ez\\, 
the system in x > becomes insulating, i.e., the resulting 
conductance vanishes. 

The reflection and transmission coefficients r±± and 
t±± are obtained by solving the continuity condition at 



with 



i-x a 



E F + aXso E F - V - aXso + £E Z 
E F - <tX so E F -V + aX so - £E Z 



(14) 



(15) 



The conductance is given by G — (e 2 /h)(2 — J2 a \ r a\ 2 )- If 
Aso -C |£f| an d £\E Z \ -C \E F — V\, the corresponding r a 
tends to zero, i.e., a perfect transmission occurs, wh ich is 
known as the Klein tunneling in monolayer graphene! 13 ! 14 ! 
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with Rashba spin-orbit interaction 

(a) E F = l.lAso (b) E F = 2A SO 
np nn np nn 
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without Rashba spin-orbit interaction 
(c) E F = l.lAso (d) E F = 2A so 
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FIG. 4. Conductance in unit of e 2 jh for the one-dimensional 
junction in the presence [(a) and (b)] (Ar, = 0.7meV) and 
absence [(c) and (d)] (Ar = 0) of the Rashba SOI. (a) and 
(c) [(b) and (d)] are the lightly (heavily) doped case as E F = 
l.lAso (E F = 2Aso)- The solid and dashed lines in (a) are 
located at £E Z = 1.5Aso and £E Z = 0.5Aso, respectively, (i)- 
(viii) are the representative points for which the conductance 
is shown in Fig. [5] These are defined as follows, (i) and (v) 

V = 4Aso (np and the double-channel regime), (ii) and (vi) 

V ~ 2Aso ( n P and the single-channel regime), (iii) and (vii) 

V = (nn and the single-channel regime), (iv) and (viii) 

V — — 2Aso (nn and the double-channel regime). 



Here we go back to Fig. [4] In the inner region of 
\Ep — V\ < ei = |Aso — £ \Ez\\, the conductance van- 
ishes since the transmitted side (x > 0) is insulating. In 
the central region of t\ < \E F — V\ < e 2 = A so + I \E Z \, 
there is a single energy band at the Fermi level, hence 
the maximum value of resultant conductance is e 2 /h. On 
the other hand, in the outer region of |£"f — V\ > 62, two 
energy bands are located at the Fermi level. Here the 
conductance becomes larger (almost double) than that 
in the central region. We refer to these regions as in- 
sulating, single-channel, and double-channel regimes, re- 
spectively. This behavior originates from a peculiarity of 
silicene, i.e., band gap and spin-split energy bands ow- 
ing to SOI and electric-field effect in the buckling struc- 
ture. Graphene, in contrast, does not have SOI nor the 
buckling structure. The resulting conductance is always 
2e 2 /h. 

For the lightly doped case [E F = l.lAso, Figs. |4^a) 
and |4jc)], one can see asymmetry of the conductance 
with respect to IE Z = Aso • To be more explicit we show 
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FIG. 5. Conductances for the one-dimensional TP/NTP 
{IE Z = 1.5A SO ) and TP/TP {IE Z = 0.5A SO ) junctions in 
the case of normal incidence. The parameters are the same 
as in Fig. [2] Cases (i)-(viii) correspond to those in Fig. |4ja) . 



the conductance as a function of V for E z = 1.5Aso 
(TP/NTP junction) and E z = 0.5A SO (TP/TP junction) 
in Fig. [5] For the double-channel regime (\E-p — V\ > e 2 ) 
[(i), (ivj, (v), and (viii)], the transmission probabilities 
of the TP/NTP and TP/TP junctions are almost the 
same. The transmission probability in the nn regime 
[(iv) and (viii)] is slightly larger than that in the np 
regime [(i) and (v)]. In contrast, for the single-channel 
regime (e x < \E F - V\ < e 2 ) [(ii), (iii), (vi), and (vii)], 
the conductances of the TP/NTP and TP/TP junctions 
are qualitatively different. Namely, the conductance for 
the TP/NTP (TP/TP) junction in the np (nn) regime 
(ii) [(vii)] takes a larger value than that in the nn (np) 
regime (iii) [(vi)]. Note that the gate- voltage dependence 
of the conductance for the TP/NTP junction is distinct 
to that for the TP/TP junction. 

This asymmetric behavior of conductance stems from 
the sublattice and spin states of the incident and trans- 
mitted electrons. Figure [6] shows the energy bands for 
x < and x > 0. The sublattice and spin states \a/3) 
for each energy band at k = are also denoted in Fig. 
[6] The incident states (x < 0) with a positive energy 

are approximately given by \-\ — ) and | ). When the 

transmitted states is given by | ) [(ii) and (vii)] or 

|H — ) [(iv) and (viii)], the conductance is large, due to 
matching of the sublattice and spin states. Especially, 

for (E F - V)/E F 0.5 and (E F - V)/E F - 0.5, the 

transmission probabilities are unity since the sublattice 
and spin states of the incident and transmitted electrons 
coincide with each other. In contrast, when the trans- 
mitted state is given by the mismatched state |++) [(iii) 
and (vi)] or | — h) [(i) and (v)], the corresponding conduc- 
tance is suppressed. Thus the matching/mismatching of 
the sublattice and spin states gives a larger/smaller con- 
ductance. 

It is emphasized that the transmitted states |++) and 
I ) are controlled by E z . As shown in Fig. [T] the two 
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FIG. 6. Energy spectra for the incident (x < 0) state with IE Z — (the left-upper and left-lower panels), and for transmitted 
(x > 0) states with £E Z = 1.5Aso [(i)~(iv)] an d with £E Z — 0.5 [(v)-(viii)]. Cases (i)-(viii) corresponds to those in Fig. [4] \otj3) 
with t z — a sgn (E z ) and a z — /3 sgn (E z ) denotes the sublattice and spin state for k = 0. 





np nn 


TP/NTP 


large small 


TP/TP 


small large 



TABLE I. Magnitudes of conductances in the silicene junc- 
tions for the single-channel regime (ei < \Ef — V\ < £2). In 
the junction, the two cases are realized: the gated region is 
nontopological (TP/NTP) and topological (TP/TP). 



states |++) and | ) are interchanged in the different 

topological phases, which is determined by E z . In other 
words, the conductance is highly tuned by E z through 
changing the symmetry of the wave function. The ob- 
tained results are summarized in Tabic. |TJ The same be- 
havior is studied in Ref. 321on the surface of a topological 
insulator with ferromagnets. 

As explained above, the conductance controlled by E z 
stems from the sublattice and spin states, i.e., the Dirac 
mass term. Therefore, this behavior does not appear in 
the heavily doped case (|-E F | > Mo,&\E z \) [Figs. |4jb) 
and Qd)] , where the mass gap ~ Aso is negligible as 
compared to E-p. The conductance asymmetry peculiar 
to the topological phase can be expected for other topo- 
logical insulators, provided that the system has a single 
Fermi surface in the pn junction. 



3. Digital field- effect transistor 

A heavily doped silicene can be applied to a new elec- 
tronic device: a silicene pn junction operates as a FET 
with digital values of conductance. From Figs. [4] (b) and 



(d), the transmission probability is almost quantized to 
for the insulating case (\Ep — V\ < ex), to 1 for the 
single-channel regime (ex < \Ep — V\ < e 2 ), and to 2 for 
the double-channel regime (\Ep — V\ > £2)- This is a 
consequence of the Klein tunneling of Dirac fermions: A 
massless Dirac fermion can tunnel through any barriers. 
And yet th e normal incident transmission probability is 
un j^yUj3LfJ] Although silicene has a finite energy gap, a 
perfect transmission approximately occurs for a heavily 
doped case, since the energy gap is effectively ignored as 
compared to the incident energy. In contrast, a graphene 
pn junction always shows a perfect transmission, i.e., the 
value of conductance is always 2e 2 /h. Graphene cannot 
be used as a FET. 



B. Two— dimensional case 

Next we turn to the two-dimensional case, which cor- 
responds to an actual silicene pn junction. The Hamilto- 
nian of the two-dimensional system in the vicinity of K 
point reads 

H(x) = -iv F d x T x - k y Ty - \ S oT z a z + £E z 0(x)t z 

- a\n(-id x a y - k y a x )r z + V6(x). (16) 

Here we assume translational invariance along the y-axis, 
i.e., the y-component of momentum k v is regarded as a 
parameter. We solve the scattering problem in the same 
way as in the previous section. 
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FIG. 7. Normalized conductance G/Go averaged over the 
incident angles for the lightly (a) (E F = l.lAso) and the 
heavily (b) (E F = 2Aso) doped cases. 



The normalized conductance G/Gq is given by 



G 
Go 



^ ~co S 9[T + (9)+T_(9)], (17) 

-tt/2 / 



with T±(8) being the transmission probability for an in- 
cident angle 9 defined by k y — fcp sin 9 and for incident 
spin ±. In the case of perfect transmission (T±(6*) = 1), 
the resulting conductance takes the value of G = 2Go, 
where factor 2 means that the system has two inci- 
dent states u + and «_ with different spin states. Here, 
Go = (k F W/7r)e 2 /h, k F = [(E F - A§o)/(«| + a 2 \l)]V*, 
and W being the width of the system. Also, Fano factor 
F, which corresponds to the shot noise-to-signal ratio, is 
given by 



F 



Gp s—y 

~G~ 



^ 2 de 

-tt/2 2 



cos 9 T a (6) [I -T a (6)]. (18) 



Figure [7] shows the normalized charge conductance 
G/Go as a function of gate voltage [(Ep — V)/E F ] and 
electric field (£E z /Xso)- Obviously, the charge conduc- 
tance of the two-dimensional pn junction (Fig. [7]) and the 
transmission probability of the one-dimensional pn junc- 
tion [Figs, [ga) andQb)] are almost the same, except 
for the broadening of line shape. This is because trans- 
port is determined basically by the normal incidence. An 
integral over the incident angle 9 solely gives line broad- 
ening of the charge conductance from that for the normal 
incidence T{9 = 0). 

The Fano factor of the two-dimensional junction is 
shown in Fig. [8] Overall, the resulting Fano factor 
for the lightly doped case (a) is smaller than that for 
the heavily doped case (b). Also, the Fano factor is 
roughly given by the inverse of conductance, i.e., it takes 
a small (large) value when the corresponding conduc- 
tance is large (small). This behavior is realized if the 
shot noise power is almost independent of the parame- 
ters (V and E z ). On the other hand, in the single-channel 
regime, the Fano factor is strongly suppressed when the 
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FIG. 8. Fano factor for the lightly (a) (E F = l.lAso) and the 
heavily (b) (E F — 2Aso) doped cases. 




FIG. 9. Silicene npn junction. The center region (the length 
L) is gated. Charge current flows along the x-axis. 



sublattice and spin states of the incident and transmit- 
ted electrons are coincides with each other (denoted by 
the arrows in Fig. [8J. Note that the Fano factor in the 
insulating region (\Ey — V\ < ex) is obtained to be unity, 
although it is not well-defined for metal-insulator junc- 
tions because the corresponding conductance vanishes 
(F —7- 0/0). We have concluded F = 1 in the insulat- 
ing region since F = 1 has been obtained for the long 
junction limit of the npn junction, as discussed in the 
next section. 



IV. SILICENE NPN JUNCTION 

Next we investigate charge transport in a silicene npn 
junction, where electrostatic field is applied in < x < L, 
which is illustrated in Fig. [9] The scattering problem of 
the npn junction is solved in a manner similar to that of 
the pn junction. The wave function ip±(x) is given by 

ip±(x = -0) = u±(kj; 0) + r± + u + (—ki; 0) 

+ r±_u_(-fo;0), (19) 



4 

£' 

i=l 



V±(0 < x < L) = ^w ±l u ai { qi ;E z )e^ x , (20) 

ip ± (x = L + 0) = t± + u + (ki; 0)e iklL 

+ t ± _M_(fc i; 0)e 4feli , (21) 
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FIG. 10. Conductance in unit of e 2 /h for the one-dimensional 
npn junction (8 = 0). L and a denote the length of the gate 
in the junction and the lattice constant, respectively. The 
incident energy is taken to be Ef = l.lAso- 



where qi and tVi are solutions of E-p — V = E ai (qi,k y ). 
Coefficients r±±, w±i, and t±± are obtained by solving 
the following boundary condition: 



V>±(-o) 

V±(i-0) 



V±(+o), 
V±(i + o). 



(22) 
(23) 



The normalized conductance and Fano factor are ob- 



tained by Eqs. (17) and (18), respectively. 



A. One— dimensional case 

First we show the conductance for the one-dimensional 
case in Fig. 10 i.e., T+(0) + T_(0). A resonant tunnel- 
ing occurs for g^L = 2mr,n G Z in a npn junction. In 
the short junction limit (L — > 0), a perfect transmission 
always occurs even when the central region is insulat- 
ing. In the short but finite- length junction [Fig. |l0[a)], 
the number of resonant peaks (tfr = 2nir/L) is still small 
(two peaks). On the other hand, the peak width is broad 
(~ 2Aso) since the length of the junction is short so that 
the transmission probability is large. Thus, two broad 
resonant peaks appear in Fig. [lO^a). As one increases 
L, the number of resonant peaks increases and the peak 
width becomes narrower, as shown in Figs. 10 b) and 
10 'c). Finally, the transmission probability of the long 
junction (L > 10000a) [Fig. [To^d)] asymptotically con- 
verges to that of the pn junction [Fig. |4][a)]. 



(a) L = 100a 
np | nn 



(b) L = 500a 
np | nn 



(c) L = 5000a 
np | nn 
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FIG. 11. Normalized Conductance G/Go in the two- 
dimensional npn junction for Ep = l.lAso- 



B. Two— dimensional case 

Next we show results on the two-dimensional case. The 
charge conductance is shown in Fi g. |11| As in the case of 
the pn junction discussed in Sec. |III B[ integral over the 
incident angle entirely causes broadening of detail struc- 
tures in the conductance. Namely, the conductance of 
the two-dimensional junction [Figs. 



11 |a)-(d)] is almost 



the same as that of the one-dimensional junction [Figs. 



10 



The Fano factor (Fig. 

TB: F 



a)-(d )]. 

In addition, wc show the Fano factor of the two- 
dimensional junction in Fig. 
[l2| is basically given by the inverse of G (Fig 
takes a small value for a resonant tunneling case. In the 
long junction limit [Figs. [l2{c) and[l2|d)], F of the npn 
junction tends to that of the pn junction [Fig. [sfa)]. And 
also, in the insulating regime {\Ep— V\ < ei), F converges 
to be unity. Namely, the Fano factor is interpreted to be 
unity for the insulating regime of pn junction. 



V. SUMMARY 

We have studied charge transport in the pn and npn 
junctions of silicene. In silicene, the topological phase 
transition occurs by applying electric field owing to the 
buckling structure. This transition affects the charge 
transport for the single-channel regime, i.e., the result- 
ing conductance is suppressed in the np regime for the 
TP/TP junction, while it is suppressed in the nn regime 
for the TP/NTP junction. We have shown that this sup- 



(a) L = 100a 
np | nn 




(c) L = 5000a 
np | nn 



(b) L = 500a 
np | nn 



(d) L = 10000a 
np nn 
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12. Fano factor in the two-dimensional npn junction for 
LlAso. 



The corresponding eigenvector |fc±) CT is given by 

\k±) s cx [g x (k) - ig v (k)} |t> + [±g(k) - g z (k)} \±). 

(A3) 

As a result, the partially diagonalized Hamiltonian 
H±(k) is given by 

H ± (k) = v F (k x T x - k y T y r) z ) + [±g(k)ri z + £E z ]t z . (A4) 

Thus, the energy spectrum is obtained to be ±E± (fc) 
with 



E 2 ± (k) =v F k 2 + [±g(k)+£E z Y 



(A5) 



The corresponding eigenvector u±(k) is given by the di- 
rect product of the eigenvectors of the pseudo spin and 
the spin as 



u ± (k) = \k±) T \k±) a 



(A6) 



with 



\k±) T oc v F (k x + ik y r] z ) |A) 

+ [E ± (fc) - (±g(k) Vz + IE Z )\ |B) . (A7) 



pression originates from matching/mismatching of the 
spin and sublattice states of the incident and transmit- 
ted electrons. Furthermore, the silicene pn junction has 
been shown to be a digital FET. It is not the case in the 
graphene pn junction, which has no band gap. The sil- 
icene junctions can be a potential new device controlled 
by two types of electric field V and E z . 
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Appendix A: Wave function in bulk silicene 

The low-energy Hamiltonian Eq. ([!]) is rewritten as 

(fc) = v F (k x r x - k y T y r\ z ) + \g(k) ■ crrj z + £E Z ] t z , 

(Al) 

with g(k) — (aX^ky, — a\^k x , — Aso)- We diagonalize 
this Hamiltonian sequentially, i.e., diagonalizing the spin 
part (a) in the first step and sublattice pseudo spin (r) 
in the next step. The eigenvalue of g(k) ■ s is obtained 
to be ±<?(fc) with 



g(k) = y/gm+9 2 y(k)+g*(k) = yja*^ k* + X 



2 

so- 



Appendix B: Symmetry 

In this Appendix, we show the symmetry of the con- 
ductance in the silicene junction. Applying 7r-rotation 
along the x-axis, one obtains 



T x a x H(k x ,k y )a x T x = H(k X7 -k y )\ E ^_ Ez 
The eigenvector u±(k) is transformed as 

TxCx^rkC&X! fcy) = Uzp{k x , ky). 

These lead to 



"(-«)(- 



-k y ,E z - 



(Bl) 



(B2) 



(B3) 



It follows that the charge conductance, which is obtained 
by the integral over k y , a, and j3, is an even function of 
E z . 

Next, we show the relation of the conductance between 
the two valleys. Applying unitary transformation t x T] x , 
the Hamiltonian is transformed as 

T x r} x H{k)T x r) x = H(k)\ E ^- Ez - (B4) 

The eigenvector is also transformed as 

T x U ± (k) = U±(k)\E z ^-E g ,Ti*->-Vz- ( B5 ) 



From Eqs. (|B3j) and (B5|, we conclude that the conduc- 
tances contributed from K and K' points are equivalent 
to each other. 



(A2) 
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